Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN I^TeX style file v2.2) 



Irradiation models for ULXs and fits to optical data 



Christopher Copperwheat^, Mark Cropper^ Roberto Soria^'^ and Kinwah Wu 

^ Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK 
^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



Received: 



o 
o 

(N 
(N 



> 

o 
o 

6 



ABSTRACT 

We have constructed a model which describes the optical emission from ultra-luminous 
X-ray sources (ULXs), and have used it to constrain the parameters of seven ULX sys- 
tems. Our model assumes a binary nature for ULXs, and accounts for optical emission 
from an X-ray irradiated companion star and accretion disk. We apply our model to 
six different ULX optical counterparts observed with HST, and one observed with the 
ESO VLT, and determine the mass, radius and age of the donor stars in these systems. 
In addition, we obtained constraints for the black hole (BH) mass in some cases. We 
use the mass accretion rate implied by the X-ray luminosity of these sources as an 
additional constraint on the donor star, by assuming the mass transfer is driven by 
the stellar nuclear evolution. We find that in general the donors are older and less 
massive than previously thought, and are consistent with being of spectral type B. 
We discuss how these results affect our understanding of the evolution and history of 
ULXs. Where we can constrain the BH masses, we find them to be consistent with stel- 
lar mass BHs or intermediate mass BHs of order ~ lOOAfQ. We make predictions for 
future observations of optical/infrared ULX counterparts, calculating binary periods 
for different BH masses in each of the seven sources. 
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1 INTRODUCTION 

Ultraluminous X-ray Sources (ULXs) are point-like, non- 
nuclear X-ray sources which have inferred X-ray luminosi- 
ties in excess of 10"^^erg s~^, assuming that the emis- 
I sion is isotropic. Some sources have X-ray luminosities 
, even greater than 10*° erg s~^, which is far in excess of 
' the Eddington limit for a stellar mass black hole (BH) 
with mass 3-20M (;). ULXs have been obse r ved in many 
nearby galaxies (Fabbiano fc Whitel l2003l : ISwartz et al.l 
I2OO4I : iLiu fc Bregma n 2005), implying that they are pop- 
ulous. To explain the large inferred isotropic luminosities 
of these sources, it has been proposed that the accreting 
object is a black hole of mass 5 — lOOOOMp- an intermedi- 
ate mass black h ole (IMBH) (jColbert fc Mushotzkvlll999l : 
iMakishima et al. I [2OOQ ). Such a mass range fits between 
those of the two well-known populations: the stellar mass 
BHs commonly found in our galaxy and the supermassive 
BHs in active galactic nuclei. The existence of IMBHs has 
important consequences for the theories of BH formation 
and evolution. It is still a mystery how supermassive BHs 
were formed. A continuous BH mass spectrum would suggest 
an evolutionary link between the three classes, thus allow- 
ing us to better understand BH formation in the younger 
universe. 

Whether IMBHs exist or not is currently under intense 



debate. The emission from ULXs may not be isotropic. 
If the emission is coUimated or relativistically beamed 
(Ki ng et al.ll200ll : iKording. Falcke fc Markogi2002l : iFabrikal 
the accretion rates of the systems do not need 
to violate the Eddington limit. However, the re is evi- 
denc e that some sources are truly ultraluminous l|Fabbiand 
|2004| ). An example is the ULXs associated with extended, 
diffuse Ha nebulae, which are more naturally explained 
by models with isotropic illumina tion of the interstellar 
medium by a cent ral X-ray source l|Pakull fc Mirionill2002l : 
iMiller et al.l |2003| ). Low- frequency quasi-periodic oscilla- 
tions (QP Os) in the X-ray emission from some sources 

20031 : 



(M82 X-1: 



Holmberg IX 



[Str ohmavcr fc Mushotzkv ^ 

X-1: iDewangan. Grifliths fc Rao ,2009 ; NGC 5408 X-1: 
IStrohmaver et alTl2007 ) suggest non-beamed emission and 
hence a truly high luminosity for those sources. Low- 
frequ ency breaks in the p ower density spectra o f NGC4559 
X-7 jCropper et al.l [ioM ) and NGC5408 X-1 jSoria et all 
l2004h are also consistent with accretors more massive than 
galactic BHs. However, population studies based on the lu- 
minosity function of ULXs suggest a more modest upper 
limit on the mass of IMBHs than previously thought. The 
break in the luminosity function at ^ 2.5 x 10"^° erg s~^ 
ISwartz et al.ll2004l : lGilfanov. Grimm fc Sunvaevll2004l ) sug- 
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gests the bulk of the IMBH population lies below ^ ISOM© , 
assuming accretion at the Eddington limit. 

Another argument proposed for IMBH accretors is 
based on the detection of a spectral co mponent interprete d 
as emission from a cool accretion disk l|Miller et al.ll20o3 ). 
However, this spectral argument has been disputed (e.g ., 
IStobbart. Roberts fc Wiim3l2006l : iGoncalves fc Soria"2006h . 
Alternative spectral models based on hig h-temperature, 
super - Eddington disks have been proposed (jStobbart et al.l 
I2OO6I : HE isawa et al ] I2OO3I I. It has also been suggested 
that super-Eddington accretion can occur in BH systems 
through a thin disk covered by a Com ptonizing corona 
l|Socrates fc Davis! I2OO5I : iDone fc Kubotal 12005 ). An accre- 
tion disc dominated by radiation pressure may exhibit 
strong density clumping. If the density inhomogeneities are 
on length scales much smaller than the disc scale height, 
such clumpy accretion discs could permit the ra diation lu- 
mino sity to exceed L^dd by factors of ~ 10 — 100 (|Begelmanl 
|2002| ). 

Since the X-ray data alone can be interpreted in dif- 
ferent ways, we have sought an alternative channel to 
X-ray investigations by which the nature of ULXs can 
be elucidated. The ULX may contain a stellar mass BH 
with super-Eddington accretion or an IMBH. The opti- 
cal/infrared (optical/IR) properties of ULXs will be strongly 
affected by the intense X-ray radiation field heating the 
donor star and accretion disk. The emission from the heated 
star may be the dominant optical component (as in the 
case of high mass X-ray binaries). Alternatively, theoret- 
ical modelling h as shown that the heated accretion disc 
might dominate (iRappaport. Podsiadlowski fc Pfahll l2005l : 



2 MODEL 

2.1 Irradiation in ULX binaries 



ICopperwheat et aLl2005|), a nalogous to low mass X-rav M- 
naries. In ICopperwheat et al. (2005 ) we considered such a 
binary model for ULX and derived a formulation which cal- 
culated the intensity and colour variations from the intrinsic 
and re-radiated emisison as a function of several parameters 
such as X-ray spectral hardness, black hole mass and incli- 
nation. The predictions of the calculations can be used to 
identify the optical counterpart, to determine the intrinsic 
properties of the component stars, and to determine the rel- 
ative contributions of the star and disk. We may also use 
the optical/IR emission as a diagnostic of the X-ray emis- 
sion and the accretion processes in ULXs/IMBHs. 

In recent times, there have been a number of opti- 
cal observations of ULXs with the Hubble Space Telescope 
(HST) and ground based 8-m class telescopes. These obser- 
vations have resulted in the isolation of optical counterparts 
of some ULXs, and have provided c olours and magnitudes 
in the optical wavebands (see, e.g. iKaaret. Ward fc ZezasI 
I2OO4I : iLiu. Bregman fc Seitze"3 l2004l : iKuntz et al.ll2005^ . In 
this paper, we apply our model to determine the properties 
of the ULX counterparts observed by these telescopes. 



In Section [2] we briefly reprise the ICopperwheat et al.l 
l|2005l) model and describe the recent refinements. In Section 
[3]we apply our model to seven ULX sources. In SectionUwe 
discuss the implication of our results on the current views 
regarding the nature of ULXs and IMBHs. 



The lCopperwheat et al" I l|2005h model assumes a binary sys- 
tem with a BH accreting matter from a companion star. We 
assume the X-ray emission from the ULX is isotropic. Our 
model considers the effects of radiative transport and radia- 
tive equilibrium in the irradiated surfaces of both the star 
and an a c cretio n disc, using the prescription described in 
IWu et al.l (|200ll ). 

We assume the disc is geometrically thin, and for the 
star we include both gravity and limb darkening effects, and 
the effects of radiation pressure. We assume the companion 
star is filling its Roche lobe and the BH is accreting matter 
via Roche lobe overfiow. This provides a geometrical con- 
straint on the system that is key in our determination of 
the optical luminosity of both the star and the disc. These 
assumptions are valid when the X-ray luminosity is in ex- 
cess of 10'"' erg s~^, since the donor star cannot generate 
a sufficient stellar wind in order to maintain the accretion 
rate required by this X-ray luminosity. However, some of the 
sources we detail in this section have X-ray luminosities of 
the order of lO'^^erg s~^. At this point it is conceivable that 
the accretion mechanism could be a wind from a massive 
supergiant. In this case we cannot constrain the binary sep- 
aration without temporal observations that reveal the binary 
period, which means we cannot quantify the stellar irradi- 
ation. However, a wind-fed BH would only intercept a few 
percent of the stellar wind from a star, and so given the mass 
accretion rate of ~ 10~^MQ/yr necessary to provide the ob- 
served X-ray luminosity, a wind loss of ~ 10~^Mq /yr would 
be required from the star, which is very high. ULXs also 
tend to have lower column densities than th e typical wind- 
fed s ources observed in the Galaxy/LMC l|Stobbart et al.l 
I2OO6I ). We therefore consider Roche lobe overflow as the ac- 
cretion mechanism to be a reasonable assumption for all of 
the sources down to lO'^ ^ erg s~ ^. 

ICopperwheat et al] l|2005l ) showed that the results of 
the irradiative calculations could depend sensitively on the 
assumed hardness of the X-rays, which is specifled by a 
parameter ^, where ^ = Sh/Ss and Sh and Ss are the 
hard and soft componen t s of t he X-ray flux respectively. 
As in ICopperwheat et al.l l|2005h , the absorption coefficients 
of the soft and hard X-rays are ksK and fc/iK, where k 
is the absorption coefficient of the optical radiation and 
ks > 1 and kh < 1 defines our soft/hard X-ray convention. 
The band boundary is a parameter to be determined. In 
ICopperwheat et al.l (|2005l ) we found that if we set kg — 2.5 
and kfi = 0.01, then for an input spectrum consisting of 
a blackbody and power-law component, this results in a 
soft/hard band boundary of 1.5keV. 

Determining the hardness ratio of the incident X-rays 
on the accretion disk is not straightforward and is compli- 
cated by the presence of absorption. Absorption tends to 
harden the X-rays. If the absorbing region which produces 
this hardened spectrum is intrinsic to the X-ray emitting re- 
gion itself, the disk and star will be irradiated by the same 
hard X-ray spectrum that we observe. If the absorbing re- 
gion is located between the binary system and the observer, 
then the irradiating X-rays will be much softer than is ob- 
served. 

By examining X-ray observations from ULX X-7 in 
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NGC 4559 (|Cropper et al.ll2004l ') we see that for this source, 
reasonable physical values for the hardness ratio range from 
5 ~ 0.1 to ~ 1. In this paper, we set the hardness ratio to 
0.1 since we expect a locally soft irradiating spectrum. In 
Section [4.61 we examine how varying this parameter affects 
our results. 

In order to determine the mass accretion rates and irra- 
diating flux for each source, we use the characteristic average 
values of X-ray luminosity as given by the authors listed in 
Table [1] We note that this may be a source of uncertainty 
in our results: since the X-ray observations are not contem- 
porary with the optical observations, the X-ray luminosity 
and spectrum at the time of the optical observation may be 
different. 

We assume the compact object in the binary system is a 
BH, but we do not assume it to be an IMBH. For all sources 
we use BH masses of 10 - IOOOMq in our model, a range 
that encompasses both a stellar mass and an intermediate 
mass nature for the BH. We do not use BH masses of beyond 
IOOOM0 since at this point the optical emission is dominated 
by a large disc, and increasing the BH mass further has little 
additional effect on our model results. Similarly, results for 
a BH mass of lOAf© can be taken to apply to BHs with 
masses less than this. 



2.2 Inclusion of stellar evolution models 



As well as detailing our model in lCopperwheat et al. I (|2005h . 
we applied it to the candidates for the optical counterpart 
of ULX X-7 in NGC 4559. We took the masses, radii and 
luminosities of various m ain sequence (MS) and supergiant 
stars from the tables in lAllenI (|l973l ). We included these 
parameters in our model, and hence were able to calculate 
the optical luminosity of the ULX if the companion star in 
the binary system was described by those parameters. We 
plotted tracks on colour magnitude diagrams, which describe 
the effect on the optical luminosity when the BH mass was 
varied. We compared these tracks to the measured colours 
and V band magnitudes of the candidates. This allowed us 
to constrain the properties of possible candidates. 

In this paper we have further developed our 
method by using t h e Gen eva stellar evolution models of 
iLeieune fc Schaered l|200lf ) to provide more realistic stel- 
lar parameters, and also to provide constraints on the de- 
rived ages. We primarily use the high mass-loss tracks which 
are recommended for us e when dealing with massive stars 
l|Maeder fc Mevrietlll994| ). These extend up to a stellar age 
of lO'^'^yr. We supplement these with the standard tracks 
when we wish to consider less massive stars with an age 
greater than this. We input the set of stellar parameters 
at each point along the evolutionary tracks into our model. 
We use our model to produce colours and magnitudes ap- 
propriate for the irradiated star and disc. We do this for 
all relevant tracks, which enables us to include stellar age 
as an additional parameter in our analysis. We repeat this 
process as we vary the other important parameters, such 
as X-ray luminosity, the inclination and orientation of the 
binary system, and BH mass. We described the parame ter 
space relatively completely in ICopperwheat et al.l l|2005h . 

An additional parameter, introduced by our use of the 
Geneva tracks, is the stellar metallicity. We use a sub- 
solar (Z— 0.2Zq) metallicity throughout. This is appropriate 



given that many ULXs are in low-metal environments such 
as dwarf g alaxies. Low-metal stars also lose less mass in stel- 
lar winds (|Eldridge fc Vinkll2006l ). Therefore, we speculate 
that they may end their lives with big ger cores, which c an 
more easily collapse directly into BHs (|Heger et al.l|20 0^). 

We note that we assume the Geneva models and metal- 
licity to be correct, so we have not taken into account any 
systematic error introduced by a difference between the stel- 
lar evolution as described by those models, and the evolution 
of a donor in a ULX binary system. A binary evolution code 
is not necessarily more appropriate, since we do not know 
when the binary evolved to a semi-detached state, and so we 
do not know to what extent the normal evolution of the star 
has been disrupted by the mass transfer. We therefore use a 
single-star model in this paper to describe the simplest case, 
where the donor star has not been significantly perturbed 
by the mass transfer. 



2.3 Mass accretion rate as a constraint of the 
system parameters 

For some of the sources in our analysis optical data alone is 
sufficient to determine the parameters of the donor star with 
good accuracy. For the others there is a large uncertainty 
even in establishing the luminosity class of the donor star, 
because the data can be fitted by a wide range of stellar 
parameters. In this situation, we can use addition conditions 
to constrain the parameter space. 

We consider a method that makes use of the information 
provided by the X-ray data and model stellar evolutionary 
calculations, which is essentially independent of the opti- 
cal/IR photometric observations. The X-ray luminosity Lx 
of an accreting object is given by the relation 



■qMc 



(1) 



where M is the mass accretion rate and rj is the efficiency 
parameter. For accretion onto BH, we may take rj — 0.1, 
which is sufficient for the purpose of this study. A ULX 
luminosity of 10*" erg s~^ therefore corresponds to a mass 
loss rate of M ~ 1.8 x 10~^MQ/yr for the donor star, if we 
assume all mass outffow from the star will be accreted onto 
the BH. 

We suppose that mass transfer occurs due to the donor 
star overfilling its Roche-lobe, and the transfer is driven by 
a gradual volume expansion of the star as it evolves away 
from its main-sequence phase. For a quasi-steady state, the 
rate of mass loss from the donor star M2 is given by 



M2 = 



Ma 



C23 — Czr 



2_ 

TJ 



_1_ 



1 



(2) 



where M2 is the mass of the donor star, C,2s and C,2r are 
the adiabatic indices of the mass donor star and its Roche 
lobe respectively, rj is the time scale of orbital angular mo- 
mentum loss, Tth is the thermal time scale of the donor star 
and Tnuc is the nuclear evolutionary time scale o f the donor 
star (|Ritteij|l98§ : lD'Antona. MazziteUi fc Ritterf l989'l. Sta- 
ble mass transfer occurs when C,2s — Car > 0. 

The orbital evolution of binary undergoi ng mass tr ans- 
fer is governed by the following equation (see|^^[l993): 



2:^-2^ 
J M2 



2\l + q) 



,(3) 
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where a is the orbital separation, J is the orbital angular 
momentum, q is the mass ratio M\ /M2 , /3 is the fraction of 
mass loss from the donor star accreted onto the BH, a is 
the specific angular momentum carried away by mass loss 
from the system, and "dot" denotes the time derivative. For 
conservative mass transfer, which is assumed in this study, 
we have J = a = 0, and /3 = 1. It follows that 



a ^ /q-l\M2 
a \ q ) M2' 



(4) 



The Roche-lobe radius Rl of the donor star and the 
orbital separation a is well aproximated via 



0.49(7 



-2/3 



(5) 



B±_ _ 

a ^ 0.6g-2/3 +ln(l + g-i/3) 

l|Eggleton|[l98^ ). The Roche-lobe filling condition requires 
Rl = R, where R is the radius of the donor star. Suppose 
that the mass transfer is quasi-steady, i.e. Rl ~ R. By com- 
bining these conditions with Equations |3] and O we obtain 



M2 
M2 



R 
R 



2(l-g) + (l + g) [2/3-5(g)l 



where the function q{q) is given by 



(2/5)g 



-2/3 



-1/3 



[3(l + g-V^)]-^ 



(3/5)g-2/3 +ln(l + g-i/3) 



For mass transfer driven by nuclear evolution, 



Ma 



-7- 



M2 



(6) 



(7) 



(8) 



where 7 is a positive, slowly varying parameter depending 
on the spectral type and mass of the donor star and the 
orbital parameters. The faster the donor star evolves, the 
higher the mass transfer rate will be, and how rapid mass 
transfer occurs is dictated by the evolutionary timescale. 
The rate of expansion of the stellar radius R and hence the 
mass accretion rate M at particular evolutionary stage of the 
donor star can be derived from stellar evolutionary tracks. 

Unless stated otherwise, we set the the mass accretion 
rate for sources with X-ray luminosities Lx ~ 10*" erg s^^ 
to be ~ Mq /yx and other sources with Lx ~ 10'^^ erg 
s"Ho be ~ lO"^M0/yr. 



2.4 Mass transfer timescale considerations 

Some of the models constructed to support stellar-mass BHs 
in ULXs propose they are intermediate- or high-mass X- 
ray binaries undergoing a phase of mass transfer on thermal 
timescales ([KinE ot al. 2001; King 2002). 

Thermal-timescale mass transfer can occur when the 
ratio q of donor mass to BH mass exceeds a critical value 
qcrit ~ 1. The mass transfer is initiated by orbital evolu- 
tion of the binary or nuclear evolution of the star, and for 
mass ratio of g > qcrit the Roche lobe of the secondary will 
shrink in response to the transfer. This will result in more 
of the donor spilling out of the Roche equipotential and a 
phase of rapid mass transfer. In response to the mass loss, 
the donor will decrease in size. In general, we find that the 
thermal timescale is much shorter than the nuclear evolu- 
tionary timescale for early type stars. On the other hand, 
relatively massive donors are required to achieve q > qcrit ■ If 



a large fraction of ULXs are found to be undergoing thermal- 
timescale mass transfer then the parameter space is much 
smaller than we examine, since the runaway scenario can 
only occur with a mass ratio q > qcrit. 

In Equation [8] we take the mass transfer to be proceed- 
ing on the nuclear evolutionary timescale; thus we do not 
consider thermal-timescale transfer in this paper explicitly. 
Our model breaks down when the thermal timescale is com- 
parable to the nuclear evolutionary timescale of the donor. 
In this case the star cannot adjust to the mass loss fast 
enough to prevent runaway, unstable mass transfer. Then 
the outer layers of the star cannot remain in radiative equi- 
librium and the assumptions of our model are violated. 

Note from Equation [2] that our formalism is versatile 
enough to consider thermal timescale mass transfer (in this 
case the mass transfer rate cannot be deduced directly from 
the stellar evolutionary track). Our irradiation model is ap- 
plicable irrespective of the mechanism and timescale of mass 
transfer. 



2.5 Inclination and orientation considerations 

For all ULXs the inclination with respect to the observer is 
currently unknown. At one extreme, the orbital plane of a 
binary is perpendicular to the plane of the sky (cos(j) = 0.0). 
In this case with a thin disk, all of the optical flux we ob- 
serve will be from the star. For any other inclination the 
optical flux will also contain a disc component, the relative 
contribution of which will increase as cos(i) is increased to 
1.0. The phase of the companion star at the time when the 
observations were made is also unknown. If the star is in 
superior conjunction with respect to the observer, the ob- 
servation will be of the irradiated hemisphere of the star. 
If the star and the observer are in inferior conjunction, the 
observation will be of the hemisphere facing away from the 
BH and therefore the flux from the star will contain little or 

no irradiated com ponent. 

We noted in ICopperwheat et ahl l|2005h that the geo- 
metric constraints of the binary system determined whether 
reprocessed light from the disc or the star dominated the op- 
tical/IR emission. The disc is truncated by tidal forces, and 
so when the separation between the star and the BH is large, 
the disc is also large and hence more likely to be the domi- 
nant optical component. A large separation is a consequence 
of assuming a high BH mass, so a high BH mass implies disc 
dominated optical/IR emission. We would expect that incli- 
nation would therefore dominate the geometrical effects in 
high BH mass systems. A low BH mass generally implies the 
emission is dominated by the donor star. We would there- 
fore expect the phase of the star to have a significant effect 
on the results only in the cases where we assume a low BH 
mass. 

We concentrate our analysis on the general case where 
there is an irradiated component from both the star and the 
disc. We therefo re start by assuming cos( i) = 0.5. As we 
demonstrated in ICopperwheat et al.l (|2005l ) this orientation 
already results in a strong contribution to the optical flux 
from the disc, and increasing the inclination to cos(i) = 1.0 
has little additional effect. By considering the cases where 
the star is in superior and inferior conjunction we examine 
the phases where the stellar contribution to the observed 
optical emission is strongest and weakest for this inclina- 
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tion. We then additionally consider superior conjunction and 
cos(i) = 0.0, describing the case when all of the optical flux is 
from the irradiated hemisphere of the star. When we assume 
inferior conjunction and take cos(i) = 0.0, we would see only 
the unirradiated hemisphere of the star. In this case it would 
be appropriate to use a standard unirradiated star. These 
extreme states allow us to explore relatively completely the 
inclination and phase orientation parameter space. 



3 APPLICATION TO OBSERVATIONS 

We now apply our calculations to the observed optical coun- 
terparts of six ULXs observed with HST, and one observed 
with the ESO VLT telescope. The photometric values we 
have used are in Table [T] Where necessary, we have con- 
verted the values to absolute magnitudes using the distances 
given in the table. 

We have corrected the values for the NGC 4559 and the 
NGC 1313 ULXs using the Galactic E(B-V) v alues given in 
ISoria et al.l (|2005l ) and lMucciareUi et all (|2005l ) respectively, 
and Av/E{B - V) = 3.1. For th e MlOl ULX we use the 
E{B-V) from lKuntz et al.1 l|2005h . which includes Galactic 
reddening and reddening from the disk of MlOl. In the case 
of the NGC 5204 ULX, we follow iLiu. Bregman fc SeitzeJ 
l|2002D for the reddening correction, using uh = 10 cm 
and assume the Galactic relati on uh ~ 5.8 x 10'^^ E{B — V) 
ijBohlin. Savage fc DraQi 19781 ). The remaining objects were 
already corrected for Galactic reddening. 

We determine the uncertainty in the values given in Ta- 
ble [T] by taking the error given for the original photometric 
values, and combine this with any additional uncertainty 
we have introduced by correcting for reddening. We base 
our confidence intervals on these errors, and fix the X-ray 
luminosity and hardness. We assume the distances given in 
Table [T] to be correct, and so do not introduce any addi- 
tional uncertainty where we have converted from apparent 
to absolute magnitudes. 

For each source, we take each evolutionary track in 
turn, and comparing the available photometric observations 
with our model calculations, calculate the variation in the 
X-squared statistic along it. We then combine our tracks, 
and determine the range of the important parameters to 
68%, 90%, 95% and 99% confidence levels. The resuhs of 
our calculations are in Figures l3l to 1 1 1 1 and Tables [2] to (5) We 
show all four confidence levels in the figures, but the values 
we quote in the tables and the text are taken at the 90% 
confidence level. 

3.1 ULX X-7 in NGC 4559 

ISoria et al.l (|2005h used HST WFPC2 observations to 
study the optical environment of ULX X-7 in NGC 4559, 
a sour ce with an average X-ray luminosity of 10^°erg 
s~^ (|Cropper et aL 2004). They found eight possible candi- 
dates for the ULX optical counterpart, listing the B, V and 
Ic standard magnitudes for each in table 2 of that paper. 

We applied our model to this source in 
ICopperwheat et al.l (|2005l ). We now do this again, us- 
ing the revised model as described in Section [5] We apply 
the exti nction correction for absorption within our galaxy 
given m ISoria et al.l l|2005h to the photometric observations. 



As in lCopperwheat et al.l l|2005l ). we find that candidates 2, 
3 and 4 are consistent with our model only when we use a 
BH mass of 2^ lOOOMp and an inclination of cos(i) = 0.0. 

ISoria et al.1 l|2005h suggested that candidate 1 was the 
most likely counterpart. The findings of ICopperwheat et al.l 
(,2005. ) are consistent with this and so for our further anal- 
ysis of this source we will assume candidate 1 is the optical 
counterpart. 

In Figures [T] and [2] we plot the confidence contours for 
the stellar age, mass and radius against the BH mass, for 
an inclination of cos(i) = 0.5 in Figure [1] and cos(i) — 0.0 in 
Figure [5] We assume the star is in superior conjunction in 
both cases. 

We see first that there is an upper bound on the BH 
mass of ~ 35M0 in the cos(i) = 0.0 case. As for the donor 
star parameters, if we examine the cos(i) — 0.5 case first we 
see that the age ranges from 10^ - lO^yr, the mass ranges 
from 5 - 2OM0 and the radius is between 9 and 30-Rq , with 
the lower radii implying a higher BH mass. For the cos(i) = 
0.0 case the age and mass range is similar. The radius lies 
between 51 and 72Rq when we assume this inclination. 

These values are consistent with candidate 1 being of 
a similar mas s and age to th e other candidates within the 
error circle (S oria et al.l l2005h. with its increased luminosity 
due to the effects of irradiative heating. 

3.2 ULX X-6 in M81 

This so urce has an average X-ra y luminosity of 2 x 10"^® 
erg s-^ l|Roberts fc Warwick! I2OOOD . iLiu et al.l l|2002l ) found 
an optical counterpart they considered unique to this ULX 
(designated NGC 3031 X-11 in that paper), and reported B, 
V and I magnitudes derived from HST ACS observations 
(Table [TJ. 

We first examine the case where this source is at su- 
perior conjunction and cos(i) = 0.5. We find that for this 
orientation our model is a poor fit to the observation at the 
90% confidence level: the irradiated disk/star are together 
too luminous to match the observation for any combination 
of star and BH. We therefore chose to adjust the hardness 
ratio of the irradiating X-ray spectrum in this case in or- 
der to fit our model to the observation. When we lower f to 
0.01, we find a good fit to the data. We show the confidence 
contours in Figure |3] We see that the stellar age ranges from 
10^ - 10^'^yr, with a lower stellar age implying a higher BH 
mass. The stellar mass ranges from 1.5 - 14Af0, and the 
stellar radius ranges from 2.5 - 8-R0. 

When we assume the star to be in inferior conjunction, 
we again require a reduced hardness ratio of ^ = 0.01 in or- 
der to obtain a good fit. We find that most of the parameter 
space that we found fitted with our model for the superior 
conjunction case is contained within that here. In addition, 
we find that for low BH masses we can fit larger stars of 
radius 10 - 17Rq with the observation. These parameters 
represent cases where the majority of the optical emission is 
from the unirradiated hemisphere of the star. The small, but 
not insignificant disc contribution, as well as the Roche lobe 
shape of the star accounts for the di fference b e tween our fit- 
ted stellar parameters and those of lLiu et al.l (|2002l ) in this 
inferior conjunction case. At higher BH masses the emission 
is disk dominated and our results are not dependent on the 
phase of the star. 
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Table 1. The photometric data we have used for the ULX in this work. Where necessary, we have converted the figures given in the 
references to absolute magnitudes and have appHed reddening corrections for galactic absorption. 



Mr. 



Mv 



Ms 



Mr 



Distance (Mpc) 



1\T r^r^ ^ecQ V 7 


1 oo _i_ n 1 o 
— / .zz ± u. iy 


— ( .yjo ± U.lu 




-D.yo it U.ID 


10 


f Soria et al. 2005) 


M81 X-6 


-4.28 ± 0.04 


-4.18 ±0.03 




-4.20 ±0.07 


3.63 


CLiu et al. 2002) 


MlOl ULX-1 


-6.19 ± 0.15 


-5.92 ±0.12 




-5.81 ±0.16 


7.2 


CKuntz et al. 2005) 


NGC 5408 ULX 


-6.4 ±0.2 


-6.4 ±0.2 




-6.1 ±0.1 


4.8 


{Subaru and HST archival data) 


Holmberg II ULX 


-6.03 ± 0.19 


-5.78 ±0.11 






3.05 


fKaaret et al. 2004) 


NGC 1313 X-2 (CI) 


-4.7 ±0.18 


-4.5 ±0.18 


-4.2 ± 0.18 




3.7 


CMucciarelli et al. 2005) 


HSTMAG 


F220W 


F435W 


F606W 


F814W 


Distance (Mpc) 


NGC 5204 ULX 


-8.51 ± 0.11 


-6.49 ±0.11 


-5.44 ±0.13 


-4.38 ±0.13 4.3 


(Liu et al. 2004) 



We examine next the case when cos(i) = 0.0 (Figure[2]). 
In this orientation, we can fit our model with the observation 
for a hardness ratio of ^ = 0.1. We note first of all that there 
is an upper bound on the BH mass of 33Mq . The stellar age 
that fits with the observation ranges from 10^'^ - 10^ '^yr, 
and the mass and radius range from 3 - 5.5 Mq and 10 - 
15Rq respectively. These values fit equally well, when we set 
(, = 0.01, since the stellar luminosity is much less sensitive 
to c hanges in t he X- ray hardness than the disk. 

iLiu et all (|2002l ) found the field stars in the vicinity of 
this ULX range in age from 1.0 x 10® - 1.0 x lO^yr. If we 
assume the donor in the ULX binary is of a similar age, we 
see from Figure[2]that the stellar parameters are very tightly 
constrained in the cos(i) = 0.0 case. In the cos(i) = 0.5 case 
(Figure [3|, we see that there is a lower limit on the BH mass 
of 20Mq if we constrain the stellar age in this way. 



3.3 ULX in NGC 5204 



HST WFPC2 and ACS observations of the optica l counter- 
part t o a ULX in NGC 5204 were described in iLiu et aD 
l|2004l . This source has an X-ray luminosity of L^, ~ 3 x 
We took the photometric measurements from 



ers_ 



iLiu et al.1 (|2004h and corrected those magnitudes for inter- 
stellar absorptio n using the Nh column density given in 
iLiu et al.1 (|2004 ). 

We find our model is a very poor fit to the observation 
when we orient the system so as to include an irradiated 
disk and/or stellar component. This poor fit is caused by 
our constraint on the mass accretion rate. When we remove 
this constraint we find our model implies a companion star 
with a mass of 60 to llOM©, a radius of 13 - 15Rq and 
an age of lO'^'^yr or less. This is a young, massive star that 
is evolving rapidly in radius and is hence transferring mass 
at a high rate, so when we set an upper limit on the mass 
accretion rate as implied from the X-ray luminosity, this 
solution is precluded. The best solution with this constraint 
applied is for a BH mass of WOOMq and a stellar age, mass 
and radius of lO'^'^yr, 52Af0 and 9Rq respectively. 

We note as an aside that, when we remove the con- 
straint on the mass accretion rate, there is an upper bound 
on the BH mass of 24OM0 when we assume an inclination 
of cos(i) — 0.5. 



3.4 MlOl ULX-1 

The source designated ULX-1 in MlOl has a peak X -ray 
luminosity of ~ 1.2 x 10^'' erg s"^ iKuntz et all l|2005h re- 
ported a unique optical counterpart and reported B,V and I 
magnitudes recorded with the HST ACS instrument. From 
these data, they deduced the companion star was a B type 
supergiant with a mass of 9 to 12Mq . They also observed no 
significant variation in the optical magnitudes over a sixty 
day period. 

When we use an inclination of cos(i) — 0.5 we find a 
stellar age of 10* - 10* '^ yr, a mass of 2 - 7Mq and a radius of 
6 - 30-Rq , with lower radius values implying a more massive 
BH. 

As in the case of the ULX in NGC 5204, we find our 
model is a very poor fit to the observation when we use an in- 
clination of cos(i) = 0.0. Again, the poor fit is caused by the 
upper bound on the mass accretion rate - when we remove 
this constraint we find a star of age lO'^'^yr to lO'^ '^yr, mass 
11 - WOMq and radius 12 - 337?© fits the observation. These 
are much looser c onstraints on the s tellar parameters than 
those reported bv lKuntz et al.l l|2005l ). By allowing emission 
from both an irradiated star and disk component, the ob- 
servation fits with a much wider range of binary systems. 



3.5 ULX in NGC 5408 

NGC 5408 contains a ULX with an X-ray luminosity of 
lO^^erg s-^ This source was initially thought to be con - 
sistent with a beamed microquasar l|Kaaret et al.1 12003| ). 
More recent studies have shown disk emission and QPOs, 
which could be interpreted as evidenc e for an IMBH 
(ISoria et al. I l2004l : IStrohmaver et al.ll2007l ). However, a soft 
disk component may also be explained by alt ernative 
scenarios that do not requir e an IMBH (Stobba rt et al.l 
l2006l : lGoncalves fc Soriall2006l ). Therefore, in this and other 
sources characterized by an X-ray soft excess, optical con- 
straints on the system parameters can help break the degen- 
eracy. 

We use the archival HST/WFPC2 and Subaru observa- 
tions to determine B, V and I photometric magnitudes for 
the optical counterpart. There are in fact a number of candi- 
dates for the optical counterpart within the Chandra error 
circle: we assume here the counterpart is the source which 
appears most luminous in the V-band HST observation. 

We examine the cos(i) — 0.5 case first (Figure [5]). We 
see that we can fit a donor star of mass 6 - 24Mq to the 



© 0000 RAS, MNRAS 000, 000-000 



Irradiation models for ULXs and fits to optical data 7 




100 

BH Mass (Mjun, 



1e+09 



1e+08 




< 1e+07 



1e+06 



100000 
1000 1 



100 

BH Mass (MjyJ 



1000 





10 



100 

BH Mass (Mgun) 



1000 



10 



100 

BH Mass (Mgu J 



1000 



Figure 1. Confidence contours for the binary parameters for the „. r. « r n- m u ^ -ii u- ■ i- i- r /■\ 

" t- j< lojure 2. As tor b igurelll but with a binary inclination oi cosli) : 



source X-7 in NGC 4559, assuming candidate 1 is the optical coun- 
terpart. We assume a binary inclination of cos{j) = 0.5, superior 
conjunction with respect to the observer and the star, and a stellar 
metallicity of Z= 0.2Zq. We use an X-ray hardness ratio of ^ = 0.1. 
The red, black, green and blue lines denote the 68%, 90%, 95% and 
99% confidence intervals respectively. 



0.0. 



observational data over the entire BH mass range. These 
stars have ages of order lO^yr and radii of 23 - 44Rq. We 
see also that when we assume a BH mass of greater than 
IOOM0, we can fit more massive (< lOOM©), younger and 
more compact stars to the observation. A very massive (> 



67 Mq) donor is also possible when we use a BH mass of less 
than 30Mq. 

In the cos(j) — 0.0 case, we see from Figure [S] that an 
upper bound of llOM© exists. The donor star has an age of 
10^ - lO'^-^yr, a mass of 6 - 15Mq and a radius of 23 - 43i?o . 
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Figure 3. Confidence contours for the binary parameters for tiie 
source X-6 in M81. We assume a binary inclination of cos(i) = 0.5, 
superior conjunction with respect to the observer and the star, and 
a stellar metallicity of Z= 0.2Zq. We use an X-ray hardness ratio 
of ^ = 0.01. The red, black, green and blue lines denote the 68%, 
90%, 95% and 99% confidence intervals respectively. 



Figure 4. As for Figure[3l but with a binary inclination of cos(i) = 
0.0 and a hardness ratio of ^ = 0.1. 



A very massive (> SOM©) donor is also possible when we 
use a low BH mass. 

3.6 ULX X-2 in NGC 1313 

iMucciarelli et al.l l|2005l ') analysed archive ESO VLT photo- 



metric data of this Lx ~ 10*° erg s~^ ULX. They found two 
possible candidates for the optical counterpart of this ULX 
within the Chandra error circle. They gave B, V and R mag- 
nitudes for the candidate designated CI in that paper, and 
V and R magnitudes for the candidate C2. This second can- 
didate was not detected in the B-band. 
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Figure 5. Confidence contours for the binary parameters for tiie „• ^ \ c t^- rn u ^ -ii u- ■ i- i- r / ■\ 

° b igure 6. As tor b igure[5] but with a binary inclination oi cos(i) : 

ULX in NGC 5408. We assume a binary inchnation of cos(i) = 0.5, 

superior conjunction with respect to the observer and the star, and 

a stellar metallicity of Z= 0.2Zq. We use an X-ray hardness ratio 

of 5 = 0.1. The red, black, green and blue lines denote the 68%, 

90%, 95% and 99% confidence intervals respectively. 



We applied our model to the photometric data of both 
candidates. We find the non-detection of C2 in the B-band 
makes it inconsistent with being an irradiated donor star 
and/or disc. We therefore take candidate CI to be the coun- 
terpart. This determination has been confirmed by a recent, 
more accurate astrometric study with the HST ACS instru- 



ment (|Liu et al.ll2007l ). We see in Figure [7] that an upper 
bound of IOOM0 exists for this source when we set the in- 
clination to be cos(i) = 0.5. If we take Figures [7] and [8] 
together we see that the donor star mass ranges from 6 to 
38Mq, with the upper bound decreasing to 23Mq in the 
cos(i) = 0.5 case. For cos(i) = 0.5 the radius is 3 - 6Rq. 
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Figure 7. Confidence contours for the binary parameters for ULX 
X-2 in NGC 1313, assuming candidate CI is the optical counterpart. 
We assume a binary incUnation of cos(j) = 0.5, superior conjunction 
with respect to the observer and the star, and a stellar metallicity 
of Z= O.2Z0. We use an X-ray hardness ratio of ^ = 0.1. The red, 
black, green and blue lines denote the 68%, 90%, 95% and 99% 
confidence intervals respectively. 



Figure 8. As for Figure[7l but with a binary inclination of cos(i) = 
0.0. 



Setting cos(i) = 0.0 results in the radius increasing to 7 - 3.7 ULX in Holmberg II 
187?0. 

This ULX has an X-ray lum inosity measured at up to 
10*° erg s"^ (|Kaaret et ai.l 120041 ') . although it is highly vari- 
able. It is associated with a diffuse, photoionised nebula, 
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Figure 9. Confidence contours for the binary parameters for tiie „■ -i a r m i i ii i ■ ■ v i- r 

" / \ J? igure 10. As for t'lgure |9j but with a binary mcunation oi 

ULX in Holmberg II. We assume a binary inclination of cos(i) = 0.5, 



superior conjunction with respect to the observer and the star, and 
a stellar metallicity of Z= 0.2Zq. We use an X-ray hardness ratio 
of 5 = 0.1. The red, black, green and blue lines denote the 68%, 
90%, 95% and 99% confidence intervals respectively. 



cos(i) = 0.0. 



whose energetics suggests that the X-ray emission from the 
accreting sou r ce is truly luminous and not strongly beamed. 
iKaaret et all (|2004l ') gives both the V magnitude of the coun- 
terpart and its (B-V) colour. 

We find the available optical data does not allow us to 
constrain the BH mass for any orientation. As regards the 



parameters of the donor star, we see in Figure [5] that if we 
set cos(i) = 0.5 the stellar age and mass are quite poorly 
constrained, with the mass ranging from 6 to 82Af0 and the 
age ranging from 10^ to lO'^'^^yr. The stellar radius is better 
defined, and fies between 4 and 127?q . 

For the cos(«) = 0.0 case (Figure I10|) . the picture is 
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more complicated. Our results suggest two possibilities for 
the donor star parameters. It can either be a star of mass 37 
- 92M0 and radius 10 - 12Rq , or a much older object with 
mass 5 - 34Mq and radius 12 - 557?0. These two possibilites 
are more tightly constrained as the BH mass increases. 



4 DISCUSSION 

4.1 Classification of the donor stars 

In Table [2] we list the Zero Age Main Sequence (ZAMS) 
stellar masses and temperatures we calculate for the donor 
stars in the seven systems we have investigated. These values 
do not account for any loss due to to the mass accretion, and 
so should be increased by an amount depending on when 
the mass transfer began. In Table |3] we list the most likely 
currently observed spectral type we obtained for each source 
. We find the donor stars are in general consistent with main 
sequence or evolved giant/supergiant stars of type B, except 
in the NGC 5204 ULX. Not suprisingly, we find the donor 
stars to be larger, less massive and older than inferred when 
irradiation is not taken into account. We can also rule out 
donors of type A or later in most cases. We now comment 
on individual sources, giving the currently observed spectral 
type in each case. 

For ULX X-7 in NGC 4559, our determination of the 
mass and radius suggest the star is a supergiant, of spectral 
type B (if we assume an inclination such that there is a 
disk component to the emission) or spectral type A (if we 
assume the plane of the disk is perpendicular to the sky - 
cos(i) = 0.0). The donor has an age and mass consistent 
with the stars in its immediate neighbourhood (jSoria et al.l 
I2OO5I). 

[lIu et all l|2002h found the photometric observations of 
the counterpart to M81 X-6 fitted with an MS 0-star. They 
noted that although the photometric data can be fit by con- 
sidering an 09 MS star the colour of the data is redder than 
would be expected. They corrected for this by assuming in- 
trinsic extinction by the dusty environment of the ULX, 
which changed their determination of the spectral type to 
an 08 star. We argue instead that the red excess could be 
explained in terms of a disc component adding to the opti- 
cal emission (as in LMC X-3. lvan Paradiis et "aLlll987^ . The 
stellar parameters we calculate instead identify the donor as 
an MS or giant evolved B-star. 

For the ULX in NGC 5408, our best fit model is a giant, 
evolved B star, but we cannot rule out the possibility of very 
massive (> 45M0) O-type donor star. When we assume the 
inclination is such that the plane of the disk is perpendicular 
to the plane of the sky, the radius implied by our model 
increases, and matches that of a B-type supergiant. 

For ULX X-2 in NGC 1 313, our calculations sug gest the 
candidate designated CI in iMucciareUi et al.l l|2005f ) is most 
likely to be the optical counterpart. The low stellar radius we 
find from the model fit suggest the donor is a main sequence 
star. The most likely spectral type is B, although a late O- 
type is also a possibility , especially when we set cos(i) — 0.0. 
iMucciarelli et all (|2005') suggest that the donor is an MS star 
of type 09 - BO, which is consistent with our findings. 

For the Holmberg II ULX, when we set cos(z) = 0.5 we 
find a tight constraint on the stellar radius but a large mass 



range. These results are consistent with the donor being ei- 
ther an MS O-star or a giant B-star. We see two distinct 
possibilities in the cos(i) = 0.0 case as well, but here the 
solutions corresponding to a lower stellar mass give a higher 
stellar radius, so the don or can be clas s ifed a s either an MS 
O-star or a B supergiant. iKaaret et al.l (|2004h suggested the 
donor was of type 04V or B3Ib, which is consistent with 
our findings. 

For the ULX in NGC 5204, iLiu et all (j2004l ) reported 
the multi-band photometry to be consistent with stars of 
type 05 V, 07 III or BO lb. They also reported HST/STIS 
far-ultraviolet spectral observations, and on the basis of 
those they suggest the star is most likely to be of type BO 
lb, although they note the spectrum does contain some pe- 
culiarities for a star of this type. The STIS data are also 
consistent with the presence of an X-ray illuminated accre- 
tion disc. Our model does not provide good fits, because we 
have applied a constraint in which the mass accretion rate 
is fixed. When we relax this constraint we find the observa- 
tion is consistent with a very massive O-type MS star, and a 
mass accretion rate of an order of magnitude or more greater 
than we expect from the observed X-ray luminosity. 

The co nclusion t h at th e donor is an MS O-star dis- 
agrees with iLiu et al.l (|2004l ). (and this would be the only 
source where we have found the donor to b e a more massive 
and compact star than originally thought). iLiu et al.l (|2004l ) 
found the photometric observations to be consistent with an 
MS or giant O-type donor star, but concluded that the star 
is a B-type supergiant is based on the additional HST STIS 
MAMA/FUV data. Specifically, the equivalent width of the 
Sim A1299 line is suggestive of a star cooler than 25000K. 
An MS O-star would be hotter than this, even on its unir- 
radiated hemisphere. If we assume our interpretation of the 
optical data is correct, we must therefore suggest that this 
line originates somewhere other than on the surface of the 
star, or that the surface layers are modified by the Roche 
lobe and material loss by accretion. 

For the ULX in MlOl, we find the donor to be a B-type 
star if the optical emission has an observable disk and a 
star component. If the disk is perpendicular to the plane of 
the sky, we have the same mass accretion constrain problem 
as for the NGC 5204 ULX. If we relax this constraint, the 
observation is consistent with the donor being a late O or B 
supergiant, or possibly an MS O star. 

4.2 Comments on mass accretion rate constraints 

In the NGC 5204 ULX we have found that the mass accre- 
tion rate as inferred from the nuclear evolution timescale 
(via the evolutionary tracks) is larger than we would infer 
from the X-ray luminosity. We can explain this in one of 
three ways. Firstly, we can assume that our interpretation 
of the optical data is correct and the donor is an MS O- 
star. This implies that the accretion efficiency 77 ~ 0.01, an 
order of magnitude less than for standard disk accretion. 
Mass outflows or advective inflows are well-known possible 
reasons for sub-standard accretion efficiency. 

The second possibility is that the BH is mostly fed by 
stellar winds rather than Roche lobe overflow. This would 
invalidate both the geometry used in our model and our cal- 
culation of the mass accretion rate, since both of them are 
dependent on Roche lobe overflow as the accretion mech- 
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anism. Given that this source is a lower- luminosity ULX 
(3 X 10"^^ erg s~^) when compared to the others in our sam- 
ple, a supergiant donor feeding the ULX via a wind is a 
reasonable possibility . Howe ver, this is inconsistent with the 
findings of iLiu et al.l (|2004| ), as they report the FUV spec- 
trum shows evidence of Roche lobe overflow. 

A third possibility is that these objects are confused 
at some level in current observations at optical wavelengths, 
with more than one star contributing to an unresolved coun- 
terpart. 

For the ULX in MlOl, we have the same problem as 
for the NGC 5204 ULX if the inclination of the system is 
set to cos(j) = 0.0. We can explain this by ruling out an 
inclination of cos(i) = 0.0 in this case. Alternatively, this can 
be exlained by either of the three possibilities detailed above. 
Given that this source has the weakest X-ray luminosity 
(1 X 10^^ erg s"^) of those in our sample, the possibility that 
this system is wind fed is relatively more likely. 

4.3 Constraining the BH mass 

A key to understanding the nature of ULXs is the determi- 
nation of the BH mass. In five out of seven systems, we can 
constrain the mass of the BH based on the optical obser- 
vations, for certain inclinations. We list these five cases in 
Table H 

If we assume the ULX X-7 in NGC 4559 has an inclina- 
tion of cos(i) = 0.0, we find an upper limit on the BH mass 
of ~ 35M0 . Analysis of the X-ra y data has suggested a lower 
limit on the BH mass of SOMrr^ l|Cropper et alj [20041 '). This 
inconsistency can be accounted for by inclining the binary 
system in our model so that the optical emission includes a 
disk component. This results in the upper limit on the BH 
mass increasing. We see in Figure [T] that for an inclination 
of cos(i) = 0.5 the upper limit has disappeared. 

When we assume an inclination of cos(j) = 0.5 for the 
ULX in NGC 1313, we find an upper limit on the BH mass 
of IOOM0 . This upper limit increases as the disk component 
is reduced, and vanishes when we set cos(i) = 0.0. This BH 
range is consistent with a BH mass of ~ lOOMp) p reviously 
inferred from the X-ray data jZampieri et al.l|2004l ). 

If we use an inclination of cos(i) = 0.0 for the ULX in 
NGC 5408, we find a maximum BH mass of llOM©. When 
we incline the system so as to include a disk component, the 
upper limit disappears. In both cases this is consistent with 
the X-ray data, which implies a BH mas s of ~ lOOM^, as - 
suming accretion at the Eddington limit (ISoria et ahlbooj ). 

We note also that the ULX in NGC 5204 has a max- 
imum BH mass of ~ 2401^0 when we assume an inclina- 
tion of cos(i) = 0.5 and a very low accretion efficiency. The 
constraint disappears when we set cos(i) = 0.0. The X-ray 
data supports the pr esence of a cool thermal disk component 
jRoberts et al.ll2005l ) which may be produced by an IMBH 
or by a stellar-mass disk cooled by other processes. Unfor- 
tunately, optical mass constraints are not strong enough to 
discriminate between the stellar-mass and IMBH scenarios. 

The fifth source for which we have a BH mass constraint 
is ULX X-6 in M81. We see that when we set cos(i) = 0.0 
there is an upper limit on the BH mass of 33M0, but if we 
set cos(i) = 0.5 we find a lower limit on the BH mass of 
2OM0 . The existence of this lower limit is dependent on the 
irradiating X-ray spectrum being softer than we have else- 



where assumed and the age of the donor being comparable 
to the of the field stars. We note also that this BH has been 
suggested to have a mass of I8M0 based on analysis of X- 
ray data \lAn et al.l I2002I '). but the authors found it to be 
model dependent. 

4.4 Predictions for further opt/IR observations 

We now look at what these new identifications for the donor 
stars mean for future optical observations of these sources. 
We show our predictions for the binary periods of the seven 
systems in Table [S] The ULX in M81 is a good prospect for 
determination of th e BH mas s, so w e base this discussion 
around this source. iLiu et all (I2OO2I ) suggested the donor 
star was of type 08V/09V. A typical 08V star has a mass 
of 23M0 and a radius of 8.57?0. Assuming a Roche-lobe 
filling donor, the binary period for such a source will be ~ 39 
hours if we assume a IOM0 BH and ~ 44 hours if we assume 
a IOOOM0BH. The amplitude of the optical lightcurve will 
be dependent on the BH mas s, binary inclination and o ther 
parameters, but in figure 3 of ICopperwheat et al.l l|2005h we 
see that the optical amplitude of the lightcurve of a typical 
O-star irradiated by an X-ray source of Lx ~ 10'*''erg s~^ is 
of the order of a tenth of a magnitude. 

As we see in Table [5] our predictions for the binary 
period of the M81 ULX change significantly with both in- 
clination and BH mass. A measured period of ~ 100 - 300 
hours would imply an inclination of close to cos(i) = 0.0 
and a low (~ IQMq) BH mass. A shorter period would im- 
ply both a more face-on system and a more massive BH, 
with a period of around 16 hours suggesting a BH mass of 
lOOOAf0 . We note also that the donor star in the cos(i) = 0.0 
/ lOM0case is a larger, more evolved object. We would ex- 
pect the variation in optical luminosity to be observable in 
this case, since the variation is not diluted by the accretion 
disk, and the optical luminosity of evolved stars is more sig- 
nifcantly modified by X-r ay irradiation than for MS stars 
jCopperwheat et al.lf2005h . A small or unobservable optical 
variation would imply an inclined system, an MS donor and 
a mor e massive BH. 

In lCopperwheat et al.l (|2005| ) we noted that infrared ob- 
servations might be an important diagnostic in determin- 
ing the nature of ULXs. We illustrate this point in Figure 
111! We show here colour-magnitude diagrams for B against 
[B — V) and H against [H — K) for the optical counterpart 
to the ULX in NGC 5408. In both cases, we plot the sec- 
tions of the evolutionary tracks which fit the optical data 
to the 90% confidence level. We plot tracks for inclinations 
of cos(i) — 0.5 and cos(i) — 0.0 with the star in superior 
conjunction. In the B versus {B — V) plot, the sets of tracks 
for the two inclinations are necessarily similar in colour and 
magnitude. However, we see that at infrared wavelengths 
there is a clear distinction between the cos(i) = 0.5 tracks 
and the cos(i) = 0.0 tracks in both colour and magnitude. 
It is clear therefore that the combination of IR and optical 
observations will allow us to constrain parameters in ULX 
systems with much higher precision. 

4.5 The evolution and history of ULXs 

If IMBHs do indeed exist, it is of great importance to clarify 
how the ULX/IMBH and the star formation process in their 
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Table 2. Calculated stellar ZAMS parameters for the donor stars. These values do not account for any loss due to to the mass accretion, 
and so should be increased by an amount depending on when the mass transfer began (* applies only when we do not apply a constraint 
on the mass accretion rate. ** applies only for a hardness ratio of g = 0.01.) 

BH Mass = 1OM0 BH Mass = IOOMq BH Mass = lOOOM© 

Mass (Mq) Log Temp (K) Mass {Mq) Log Temp (K) Mass {Mq) Log Temp (K) 



NGC 4559 X-7 



cos(i) 


= 0.5, superior conjunction 


5.08 - 


12.19 


4, 


,30 


- 4, 


,48 


cos(j) 


= 0.5, inferior conjunction 


O.Do — 


iZ.Oo 


A 




A 


A O 


cos(j) 


= 0.0, superior conjunction 


5.97 - 


12.77 


4, 


,33 


- 4, 


,49 


M81 X-6 














cos(i) 


= 0.5, superior conjunction ** 


- 




- 








cos(i) 


= 0.5, inferior conjunction ** 


- 




- 








cos(i) 


= 0.0, superior conjunction 


2.48 - 


5.37 


4, 


,11 


- 4, 


,31 


NGC 


5204 ULX 














cos(i) 


= 0.5, superior conjunction * 


68.4 - 


109.3 


4, 


,71 


- 4, 


,74 


COs(2) 


= 0.5, inferior conjunction * 


69.4 - 


110.3 


4. 


,71 


- 4, 


,74 


cos(j) 


= 0.0, superior conjunction * 


65.1 - 


116.3 


4. 


,70 


- 4, 


,75 


MlOl 


ULX-1 














cos(i) 


= 0.5, superior conjunction 


2.46 - 


4.13 


4, 


11 


- 4, 


,24 






81.3 - 


84.3 


4. 


72 






cos(i) 


= 0.5, inferior conjunction 


3.24 - 


4.03 


4, 


,18 


- 4, 


,24 






81.3 - 


84.3 


4, 


,72 






cos(i) 


= 0.0, superior conjunction * 


15.7 - 


101.3 


4, 


,52 


- 4, 


,74 


NGC 


5408 ULX 














cos(j) 


= 0.5, superior conjunction 


5.08 - 


23.49 


4, 


,30 


- 4, 


,59 






68.4 - 


117.3 


4, 


,71 


- 4, 


,75 


cos(j) 


= 0.5, inferior conjunction 


5.37 - 


14.3 


4, 


,31 


- 4, 


,51 






70.4 - 


117.3 


4, 


,71 


- 4, 


,75 


cos(j) 


= 0.0, superior conjunction 


5.08 - 


14.35 


4, 


,30 


- 4, 


,51 






79.2 - 


117.3 


4, 


,72 


- 4, 


,75 


NGC 


1313 X-2 














cos(j) 


= 0.5, superior conjunction 


9.29 - 


23.1 


4, 


,42 


- 4, 


,59 


cos(i) 


= 0.5, inferior conjunction 


6.71 - 


27.11 


4, 


,36 


- 4, 


,61 


cos(i) 


= 0.0, superior conjunction 


5.68 - 


26.11 


4, 


,32 


- 4, 


,61 


Holmberg II ULX 














cos(«) 


= 0.5, superior conjunction 


5.68 - 


81.4 


4, 


,32 


- 4, 


,72 


cos(i) 


= 0.5, inferior conjunction 


5.68 - 


88.3 


4, 


,32 


- 4, 


,73 


cos(j) 


= 0.0, superior conjunction 


5.08 - 


90.3 


4, 


,30 


- 4, 


,73 



5.08 - 


9.80 


4.30 


- 4.44 


5.67 - 


20.73 


4, 


,32 


- 4.57 


5.68 - 


9.11 


4.32 


- 4.42 


5.67 - 


20.73 


4, 


,32 


- 4.57 


3.38 - 


7.10 


4.20 


- 4.37 


6.36 - 


9.10 


4, 


,35 


- 4.42 


2.56 - 

- 


8.03 


4.12 

- 


- 4.39 


6.56 - 

- 


9.10 


4, 

- 


,35 


- 4.42 


72.0 - 


95.3 


4.71 


- 4.73 












72.0 - 


95.3 


4.71 


- 4.73 


- 




- 






69.6 - 


117.3 


4.71 


- 4.75 


73.4 - 


118.3 


4, 


,71 


- 4.75 


2.48 - 


4.81 


4.11 


- 4.28 


2.93 - 


6.35 


4, 


,16 


- 4.35 


2.94 - 


4.58 


4.16 


- 4.27 


2.93 - 


6.35 


4, 


,16 


- 4.35 










45.0 - 


52.4 


4, 


,67 


- 4.69 


19.0 - 


100.4 


4.56 


- 4.74 


15.7- 


104.3 


4, 


,52 


- 4.74 


5.68 - 


18.98 


4.32 


- 4.56 


8.32 - 


57.4 


4, 


,40 


- 4.69 


84.3 - 


88.3 


4.72 


- 4.73 












5.67 - 


19.0 


4.32 


- 4.56 


8.32 - 


57.4 


4, 


,40 


- 4.69 


86.3 - 


88.3 


4.72 


- 4.73 












9.80 




4.44 




- 




- 






10.3 - 


12.1 


4.45 


- 4.48 












11.0 - 


15.3 


4.46 


- 4.52 












5.37 - 


35.5 


4.32 


- 4.65 


7.56 - 


37.5 


4, 


,38 


- 4.65 


6.71 - 


68.4 


4.36 


- 4.71 


10.82 


- 22.1 


4, 


,45 


- 4.58 


5.67 - 


68.4 


4.32 


- 4.71 


10.82 


- 22.1 


4, 


,45 


- 4.58 


4.81 - 


9.11 


4.28 


- 4.42 


10.28 




4, 


,45 




57.2 - 


88.3 


4.70 


- 4.73 


72.3 - 


84.3 


4, 


,71 


- 4.72 



Table 3. Spectral type of the donor stars. We compare the classification given by previous authors (as listed in table [TJ with our 
determination of the current and ZAMS spectral types. We give the most likely type in each c ase - we elaborate o n this in t h e text . 
(* applies only w h en we do not apply a constrai nt on th e mass accr etion rate .). References: ISoria et al. I ll2005h . i Liu et al. I 1I2OO2II . 
^Liu et al.l l|2004h . ^Kuntz et all l|2005f l. ^Mucciar elli et al.l |2b05), '^K aaret et al.l 12004:). 

Previous Our determination 

spectral type Current ZAMS 



NGC 4559 X-7 O / B supergiant^ Late B - A giant/supergiant B5 - BO 



M81 X-6 
NGC 5204 ULX 
MlOl ULX-1 
NGC 5408 ULX 
NGC 1313 X-2 
Holmberg II ULX 



08V / 09V2 
BOIb^* 

B supergiant* 

BO - 09 MS^ 
04V / B3Ib6 



B MS/giant 
O MS * 
A/B MS/giant 

B giant/supergiant (or O MS) 
B MS 

B giant/supergiant (or O MS) 



B9 - B2 
05 or earlier * 
AO - 04 
B4 or earlier 
B6 - 06 
B5 or earlier 
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-0.4 -0.3 -0.2 -0.1 0.1 0.2 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 

(B-V) (H-K) 

Figure 11. Colour-magnitude diagrams showing tlie 90% confidence contours for the optical counterpart of the ULX in NGC 5408, for 
inclinations of cos(j) = 0.5 (black) and cos(i) = 0.0 (red) with the star in superior conjunction. 



Table 4. BH mass constraints, and the inclination for which they 
apply (* applies only when we do not apply a constraint on the 
mass accretion rate. ** applies only for a hardness ratio of ^ = 
0.01.) 



NGC 4559 X-7 
M81 X-6 

NGC 5204 ULX 
NGC 5408 ULX 
NGC 1313 X-1 



< 35MQfor cos(j) = 0.0 

> 20MQfor cos(i) = 0.5 ** 

< 33MQfor cos(i) = 0.0 

< 24OM0for cos(i) = 0.5 * 

< llOM0for cos(i) = 0.0 

< lOOM0for cos(i) = 0.5 



Table 5. Predictions of the orbital periods in these systems, in 
hours, based on our determinations of the donor star parameters 
(* applies only when we do not apply a constraint on the mass 
accretion rate. ** applies only for a hardness ratio of ^ = 0.01.) 



BH mass 



lOAf, 







lOOA/0 



IOOOMq 



cos(j) = 0.5, superior conjunction 



NGC 4559 X-7 


418 - 


732 


260 - 


- 276 


67.0 - 


- 78.8 


M81 X-6 ** 






47.9 


- 59.5 


15.9 - 


- 16.2 


NGC 5204 ULX * 


27.4 


- 42.0 


37.8 


- 44.5 






MlOl ULX-1 


767 - 


873 


271 - 


- 295 


77.6 - 


- 82.9 




26.6 


- 26.9 










NGC 5408 ULX 


114 - 


277 


94.9 


- 112 


33.4 - 


- 36.6 




27.5 


- 39.3 


37.1 


- 37.4 






Holmberg II ULX 


27.2 


- 152 


35.4 


- 61.5 


21.1 - 


- 22.1 


NGC 1313 X-2 CI 


22.6 


- 36.4 


21.8 


- 27.4 







cos(j) = 0.0, superior conjunction 



NGC 4559 X-7 


864 - 


2060 










M81 X-6 


136 - 


287 










NGC 5204 ULX * 


27.5 


- 46.7 


38.8 


- 56.4 


41.8 - 


- 57.7 


MlOl ULX-1 * 


27.2 


- 387 


38.2 


- 394 


40.6 - 


- 376 


NGC 5408 ULX 


247 - 


732 


760 










27.5 


- 41.0 










Holmberg II ULX 


27.0 


- 417 


542 - 


- 1500 


783 










35.3 


- 37.3 


36.2 - 


- 38.5 


NGC 1313 X-2 CI 


26.0 


- 96.2 


28.4 


- 157.5 


28.3 - 


- 190.4 



vicinity are related. The open question to resolve is whether 
the donor star is coeval to the BH progenitor or captured by 
the BH some time after formation. If the star and the BH 
formed together, determining the age of the donor star also 
determines the age of the BH. If the star was captured by 
the BH, then the statistics of the spectral type and mass dis- 
tribution of the donor stars can be used to set constraints on 
the capture rate and hence provide estimates to the IMBH 
populations. 

In our study, we found that the donor stars are mainly 
of spectral type B, and are sign ificantly o l der th an previ- 
ously determined. For example, iLiu et al.l \2m± inferred 
the donor star in M81 X-6 ULX to have an age of less 
than lO^'^'yr. We find the minimum stellar age to be an 
order of magnitude greater, if the BH mass is assumed to 
be <~ lOOA/0. In a number of cases, the photometric data 
alone allows us to infer the donor is of spectral type B. In 
other cases, a range of (more massive) possibilities exist, but 
by applying constraints on the mass accretion rate we find a 
B-type star to be the most likely donor. In the case where the 
optical data points to a massive, O-type donor, the implied 
mass accretion rate is inconsistent with that which would 
be expected from the X-ray observation, given our assumed 
accretion efficiency of »7 = 0.1. We suggest therefore, that 
donor stars of a narrow spectral and mass range are nec- 
essary to produce a very luminous, Roche lobe fed ULX, 
and our finding of large fraction of B-type stars in the ULX 
sample may be significant. 

If the compact objects in these systems are indeed 
IMBHs, and if we assume the capture scenario, then the 
fact that a B-type donor is sufficient to fuel a ULX allows 
a lower spatial density for IMBHs for the observed pop- 
ulation of ULXs than if the donors were found to be of 
type O, since B-type stars are more common and so the 
chances of forming a ULX binary are higher. Various au- 
thors have mo delled the tidal capture of a donor star by an 
IMBH (see e.qlHopman. Portegies Zwart fc Alexandeii2004l : 
iBlecha et al.|[2006l ). The capture rate appears to be ~ stel- 
lar number density, but is only weakly dependent on stellar 
mass. This would imply that more B stars than O stars 
should be captured. However, a competing effect is the re- 
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quirement that the star comes close enough to be captured 
but far enough not to be tidaUy destroyed. This requirement 
may favour O stars, since they may more easily survive tidal 
squeezing. We note also that tidal capture of isolated stars 
is only one process through which an IMBH might aquire a 
companion; another process is by capturing stars in binary 
systems, which may have a different frequency of occurence 
and period distribution for B or O stars. 

It is also interesting to note that the two lowest lumi- 
nosity ULXs in our sample, those in M81 and MlOl, are 
also those where an old, less massive donor of age ~ lO^Myr 
is a possibility. This may be related to the existence of a 
population of low- luminosity sources (<~ 2 x 10^^ erg s~^) 
also found in old elliptical galaxies, probably identified as 
low-mass XRBs. Conversely, the ULXs more luminous than 
2 X lO'^^erg s~^ are almost always found in star-forming 
galaxies, and we find the donor stars to be of age ~ lO^Myr 
or less. 



4.6 Other systematic effects 

We have made a number of assumptions in this work. We 
reexamine here some of the effects that will have an influence 
on our conclusions. 



4-6.1 System geometry 

We have assumed a thin disk in our model, but as we noted 
in Section [1] some authors have suggested more complicated 
disk mo dels, with a thin disk covered by a Comptonised 
corona (ISocrates fc David l2005h . These models have been 
developed to support the possibility of super-Eddington ac- 
cretion in ULX systems. The corona emits hard X-rays 
which are reflected by the ionized surface of the inner disk. 
As in our disk model, the bulk of the optical emission comes 
from the outer regions of the disk, due to these regions hav- 
ing a much larger surface area than the inner parts of the 
disk. To a large extent therefore the optical/IR emission is 
similar. Ifowever, the corona model will result in a harder 
X-ray spectrum incident on the outer disk. We discuss the 
effect of changing the X-ray hardness in more detail in Sec- 
tion U^HISl Additionally, since the corona extends away from 
the disk surface the angles of incidence of the X-rays on the 
outer disk regions are changed somewhat, but the extent of 
the corona is not large enough for this to have a significant 
effect. 



4-6.2 Radiation pressure 

In the results we have presented here, we have not included 
the effects of radiation pressure. The effects of radiation 
pressure on the shape of a Roche lobe filling star are un- 
clear: some authors suggest that the shape will be unaffected 
jHow arth 1997|). whereas others predict a significant effect 
in ve ry X-ray lumin ous binaries (IPhillips fc Podsiadlowskil 
|2002|). We showed inlCopperwheat et al.1 l|2005l ') that under 
the iPhflhps fc Podsiadlowskil l|2002l ) formulation, radiation 
pressure has an increasing effect on the shape and luminos- 
ity of the donor as the BH mass is decreased. This is because 
as we decrease the BH mass the binary separation also de- 
creases, and the X-ray flux incident on the stellar surface 



increases. We found also that the effect decreases as we use 
a donor of later spectral type, for the same reasons. We find 
therefore that for most of our fitted solutions, where the 
BH mass lies between 100 and IOOOM0 and the donor is 
found to be a star of type B or later, the effect of including 
radiation pressure on our results is small. We begin to see 
an appreciable deviation from the results we have reported 
when we use a BH mass of ~ IOOMq and a main sequence, 
O-type donor. Even in this case, the fitted stellar mass only 
changes by 1 or 2Mq, which is a small percentage of the 
total stellar mass and not enough to alter our classification 
of the star. We would expect a mo re significant deviation for 
low (~ IOMq) BH masses, but the lPhillips fc Podsiadlowskil 
(2002) formulation we use for radiation pressure becomes 
inappropriate at this point, since it does not allow for any 
circulatory currents in the stellar surface driven by the ir- 
radiation. This means it tends to represent an extreme case 
of maximum stellar distortion, and in this extreme case the 
radiation pressure is very large, to the point of stripping the 
outer layers away from the donor so that it looks quite unlike 
an ordinary star. A radiation pressure formulation which in- 
cluded circulatory currents would predict a lessened effect 
on the stellar shape. Shielding by the disk will mitigate this 
effect further. 

Given that the actual effect of radiation pressure on 
the stellar shape in ULX systems is not fully understood, 
we consider our omission of this component in this paper to 
be appropriate. We do however note this introduces an addi- 
tional source of uncertainty, particularly for low BH masses. 

4-6-3 X-ray hardness 

The hardness of the X-ray spectrum determines the depth 
at which the incident radiation deposits most of the energy. 
Soft X-rays are easily absorbed at the disk surface, while 
hard X-rays attenuate only at large optical depths. For inci- 
dent X-rays with a soft spectrum, a hot surface skin layer is 
formed at depths optically thin to optical radiation, and the 
emission from the skin layer is at wavelengths shortward 
of the optical bands. However, for incident X-rays with a 
hard spectrum, most of the energy is deposited at depths 
optically thick to optical radiation. This heats the internal 
regions of the disk plane and hence leads to a more lumi- 
nous disk at the optical wavelen gths. This was illustrated 
in figure 6 of ICopperwheat et al.l (.2005), where we see that 
as the X-ray hardness is increased, the disk V magnitude 
similarly increases. We noted also that the effect of varying 
the hardness ratio on the stellar luminosity was very small. 

Because we do not know the hardness of the X-ray spec- 
trum incident on the irradiated surfaces in these systems (as 
discussed in Section [5TT| we have used a fixed hardness ratio 
of 5 = 0.1. We now discuss the effect of varying the hardness 
ratio on our results. 

In Figure [12] we have plotted our determinations of the 
masses and radii of the donor stars in the NGC 4559 and 
NGC 5408 ULXs. We have fixed the inclination to cos(i) = 
0.5, the BH mass to lOOM© and the stellar age to lOMyr. 
If we examine the results for the ULX in NGC 4559, we see 
that varying ^ has little to no effect on our determination 
of the mass of the donor star. Our determination of the 
radius however, decreases with increasing hardness ratio. We 
see a similar decrease in stellar radius with increasing ^ in 
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Figure 12. The change in stellar mass/radius with hardness ratio ^ for ULX X-7 in NGC 4559 (left) and the ULX in NGC 5408 (right). 
We use an inclination of cos(j) = 0.5, a BH mass of IOOA/q and a stellar age of lOMyr. The lines show the range of stellar parameters 
that fitted with our model at the 90% confidence level. For clarity, we mark with an 'x' cases where our model produces a single unique 
solution. 



the NGC 5408 case. We see also an appreciable decrease 
in stellar mass in this case, but this decrease is still much 
smaller than the decrease in stellar radius. 

In a Roche lobe geometry the shapes of the lobes is 
determined by the mass ratio of the two components, and the 
scale of the system is set by their separation. In our model 
we determine the scale by setting the volume radius of the 
secondary Roche lobe to be equal to the undistorted radius 
of the donor star. The results of Figure [T^] can therefore 
be understood as follows. A harder incident X-ray spectrum 
leads to an accretion disk which is more luminous at optical 
wavelengths. If the disk in our model is hotter, then to keep 
it consistent with the observation it must be smaller, so the 
scale of the Roche lobes must be decreased. Our calculations 
therefore result in a smaller fitted stellar radius. Since the 
mass ratio affects only the shape of the Roche lobes, the size 
of the disk is only weakly dependent on stellar mass, and so 
varying the hardness ratio will tend to have little effect on 
this parameter. 

In both ULXs, the changes in stellar parameters are 
smaller than might be expected, given the significant varia- 
ti on in disk luminosity with hardness ratio shown in figure 6 
of lCopperwheat et al. I l|2005l ). This is because that figure did 
not include the stellar component. Changing the mass ratio 
results in a different determination of the stellar radius, but 
changing this parameter results in a change in the luminos- 
ity of both the star and the disk. The radius therefore does 
not need to be changed by much to have a large effect on 
the overall luminosity. 

Through examination of the X-ray data, we have found 
^ to vary from ~ 0.1 to ~ 1 in the NGC 4559 ULX. We have 
assumed this to be a physically appropriate hardness ratio 
range for all the ULXs in our sample, and have used ^ = 0.1 
for all the model fits in this paper. Figure [12] indicates that 
the general findings in this paper will not be invalidated 
if the X-ray spectrum in these systems is harder than we 
have assumed. We see that increasing the X-ray hardness to 
^ = 1 causes no change in the stellar mass in one case and 
a decrease of ~ IMq in the other. We see also a decrease 
in the stellar radius of 2 - 5i?0 as ^ is increased from 0.1 



to 1. We conclude therefore that, if the X-ray spectrum is 
harder, our determinations of the spectral type of the donors 
are still valid but they may be somewhat smaller and less 
evolved than our results suggest. 

4-6.4 Perturbation of the donor as a result of mass 
transfer 

A key distinguishing feature between the IMBH evolution- 
ary scenarios is the epoch at which the steady mass transfer 
began. There are two possibilities: (a) we observe the sys- 
tem where mass transfer has not yet significantly affected the 
state of the star, and so the models we use are a fair descrip- 
tion of the donor; (b) we observe it after the mass transfer 
has already significantly altered the state of the star and a 
binary evolution code is required for proper modelling. 

We note that even in cases when a binary evolution code 
should be used the results we present here are not initially 
invalidated. Our method is to determine the star which, in 
the presence of the intense X-ray radiation field produced by 
the ULX, will have the appearance of the observed donor, 
as it exists at the current epoch. Therefore, assuming our 
model assumptions (such as the assumption that the star 
is in thermal equilibrium) are correct, the masses and radii 
will be accurate, since these parameters depend on the cur- 
rent physical state of the star, not its history. Conversely, 
the parameters we have determined such as the stellar age 
and ZAMS mass will be inaccurate if the star has lost a 
significant amount of its mass through Roche lobe overflow. 

A binary stellar evolutionary code is necessary when 
studying a ULX in which the mass transfer has endured 
long enough for the appearance of the donor to be signifi- 
cantly modified. It is unclear as to how long it takes for a 
star to deviate in characteristics from the single star evolu- 
tionary tracks when it is undergoing mass loss at the rates 
we have inferred. The mass loss can be assumed to have lit- 
tle affect on very massive stars (~ 40Mq and above) since 
over the course of their short lives they will transfer only a 
small percentage of their mass onto the BH through Roche 
lobe overfiow. On the other hand, it has been shown that 
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low mass stars (~ IM©) deviate from single star models 
very rapidly when un dergoing phases of high mass transfer 
ijSchenker et al.ll2002[ l. 

It is reasonable to expect that a star that has trans- 
ferred, say, 50% of its mass onto the BH through Roche 
lobe overflow will be significantly perturbed. When we ex- 
amine the accumulated mass loss for star with ZAMS mass 
10 - 2OM0, we find that to lose 50% takes ~ 2 - 3Myr 
when an X-ray luminosity of 10*° erg s~^ is assumed, al- 
though this varies depending on the point in the stellar evo- 
lution at which the mass transfer begins. This is a short 
length of time, but our model predicts the length of time 
in which any star can sustain mass transfer at ULX rates is 
also short, even when the efi^ect of the mass loss on the star 
is not considered. Massive stars can transfer mass at ULX 
rates from ZAMS, but have intrinsically short (< lOMyr) 
lifespans. Lower mass donors only begin to transfer mass at 
the required rate towards the end of their time on the MS. 
The maximum duration of the binary as a ULX, assuming a 
10*° erg s~^ X-ray luminosity and a donor ZAMS mass of 10 
- 2OM0 , is around 6 - 8Myr. Assuming a lower mass transfer 
rate, as would be expected in ULXs with X-ray luminosities 
of ~ 10^® erg s~^, means that the donor takes significantly 
longer to be perturbed by the mass transfer. We note again 
that these conclusions are also dependent on the accretion 
efficiency we have assumed in ULXs being correct. 

Note that in this estimation we have assumed the donor 
is unperturbed at the beginning of the ultraluminous X-ray 
emitting phase of the binary. Our model suggests that low 
mass stars can transfer mass at sub-ULX rates in the early 
part of their lives, assuming they evolve to a semi-detached 
state shortly after ZAMS. This would mean that the star 
would already be appreciably perturbed before it began its 
ULX phase. Similarly, if the star has at some point in its 
past undergone a phase of extreme mass transfer (such as 
thermal-timescale mass transfer, as discussed in Section ETI)) . 
then a binary evolution code would almost certainly be nec- 
essary. We note also that very high rates of mass trans- 
fer are possible even when the driver is nuclear evolution: 
very massive stars can transfer mass to the BH at a rate 
or more, particularly towards the end of 
their MS life and beyond. This would result in rapid and sig- 
nificant deviation from the single star tracks we have used. 
Mass transfer at this rate is in excess of what we would 
expect for the ULXs in our sample, but some very X-ray 
luminous systems such as the Lx =~ 10*^ erg s~^ ULX in 
M82 could have very early type donors in states of extreme 
mass loss. 

We conclude therefore that a binary evolution code may 
be necessary in some of the sources discussed here, but 
given that the timescales of the ULX active phase and the 
timescale for the star to be perturbed by the mass loss are 
similar, the use of single star tracks is reasonable for this 
investigation. 



5 CONCLUSIONS 

We have constructed a model which describes the opti- 
cal emission from ultra-luminous X-ray sources (ULXs). 
Our model assumes an X-ray binary nature for ULXs, 
with the donor star overfilling its Roche lobe. We as- 



sume radiative equilibrium and use a radiative transfer 
formulation to determine the re-radiated thermal emission 
ICopperwheat et al.ll2005[ ). We predict the optical luminos- 
ity of the emission from both the irradiated star and a thin 
accretion disc. We use the mass accretion rate as implied by 
the observed X-ray luminosity as an additional constraint 
on the donor star parameters, assuming that the nuclear 
evolution of the star is the driver for the mass transfer, and 
hence calculating the mass accretion rate from the rate of 
increase in radius of the star, compared to the evolution of 
the Roche lobe. 

We have applied our model to optical observations of 
seven different ULX optical counterparts, and determined 
the parameters of the donor stars, by fitting the multi-band 
photometric observations to the optical emission predicted 
by our model for different sets of stellar parameters. The 
stellar m odels were taken from the G eneva stellar evolution 
tracks of iLeieune fc Schaererl l|200ll ). We varied the other 
parameters in our model, such as binary inclination, position 
of the star with respect to the observer when the observation 
was made, BH mass and stellar metallicity, in order to better 
determine the nature of the donor star. 

Previous analysis found that the donor star of ULX 
are often MS O stars or early-type supergiants. Here we 
have shown that the donor stars are older and less massive. 
Their spectral types are generally consistent with MS stars 
or evolved giants/supergiants of spectral type B or later. 
They tend to have ages ~ 10 - lOOMyr. However the accre- 
tion rate necessary to fuel an X-ray luminosity of 10*° erg 
s~^ will result in a B-star being completely consumed in 5 
- lOMyr. This implies an upper limit on the duration of the 
mass transfer. It may be that the capture of the donor by 
the BH occured within this period of time, or it could be 
that binary formed a long time ago but evolved to a semi- 
detached state within this time. A further possibility is that 
the donor was originally a much earlier-type, more massive 
star that has been transferring mass over most of its lifetime, 
such that it has the appearance of an older and less massive 
star in the current epoch. The use of binary stellar evolution 
models will be required to fully explore this possibility. 

The ULX in NGC 5204 is an exception to our general 
findings. Our model cannot provide a good fit with the ob- 
servation in this case, owing to the constraint we apply on 
the mass accretion rate. When we remove this constraint 
the best solution is for an MS 0-star in NGC 5204, a more 
massive donor than has been suggested for this system by 
other authors. However, our fits suggest a much higher mass 
accretion rate than would be expected from the X-ray lumi- 
nosity. A possibility is that this ULX is wind fed, instead of 
via Roche lobe overfiow as we have assumed. 

There are a number of systematic effects, which we 
have discussed. It should be noted that in a number of 
cases these effects may be greater when we assume a stel- 
lar mass BH rather than an IMBH, due to the physics of 
super-Eddington accretion being uncertain. In stellar mass 
BH systems the disk geometry may be significantly different 
from the thin disk approximation we have used, and if there 
is strong beaming mass transfer may be occuring on thermal 
timescales. In addition, the effect of radiation pressure may 
significantly perturb the star, whereas we have shown this 
effect is small when a more massive BH is used. 

We have also discussed the implications for the binary 



© 0000 RAS, MNRAS 000, 000-000 



Irradiation models for ULXs and fits to optical data 19 



periods and optical amplitudes, which will be revealed by 
future observation of those sources. Our model allows us 
to determine the relative contribution of the star and the 
disc to the emission. For a less massive BH, the optical/IR 
emission will be dominated by light from the star. For a more 
massive BH, the optical/IR emission will be disc dominated. 
If we assume an irradiating X-ray luminosity of 10*'' erg s~^, 
the disc will dominate the emission for a BH mass of > 
90Mq if we assume the donor is a typical B-type MS star, 
and for a BH mass of > 3OOM0 if we assume a B supergiant. 

In more than half of these systems, we can constrain the 
mass of the BH based on the optical observations, although 
the existence of these constraints are dependent on an as- 
sumed inclination. If we assume an inclination of cos(i) = 0.5 
for the ULX X-1 in NGC 1313, we find the BH has a max- 
imum mass of ~ IOOMq. If we assume an inclination of 
cos(i) = 0.0 for the ULXs in NGC 5408 and NGC 4559, 
we find the BHs in these systems have a maximum mass of 
~ llOM© and ~ 35Mo respectively. For ULX X-6 in M81, 
we found a minimum BH mass of 20Mq, or a maximum 
mass of 33Af0, depending on the inclination. A temporal 
optical study of this source could determine with some ac- 
curacy the spectral type of the donor, the inclination and 
the BH mass, so we conclude that this source in particular 
is an excellent target for further observations with the aim 
of positively discriminating between stellar mass BH and 
IMBH scenarios. 
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